Influence of coaxial cylinders on vortex breakdown in a closed flow. 
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The effect of fixed cylindrical rods located at the centerline axis on vortex breakdown (VB) is 

qq". studied experimentally and numerically. We find that the VB is enhanced for very small values 

' of the rod radius d, while it is suppressed for values of d beyond a critical value. In order to 

I characterize this effect, the critical Reynolds number for the appearance of vortex breakdown as a 

, function of the radius of the fixed rods and the different aspect ratios was accurately determined, 

^ ^ ' using digital particle image velocimetry. The numerical and experimental results are compared 

' showing an excellent agreement. In addition, a simple model in order to show that this effect also 

, appears in open pipe flows is presented in the appendix. 

I ■ PACS numbers: 47.32.-y, 47.32. cd, 47.32. Ef 

Keywords: vortex breakdown, recirculation flow, control 
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^ , The development of structural changes in vortical flows and, particularly, vortex breakdown [l|, [3, H, 0, has been 
intensively investigated during the last years [Gj, iTj, iJ, H, [l^, [ll|, [12, 13, l4, [H, HI, 113, EBl- The characteristic and 
\ fundamental signature of vortex breakdown is the appearance of a stagnation point followed by regions of reversed 
. ^ ' axial flows with a bubble structure when the swirl is sufficiently large. This structural change is also accompanied 
^1 by a sudden change of the size of the core and the appearance of disturbances downstream the enlargement of 
^ ^7 *. the core. Vortex breakdown is very important in several applications of Fluid Mechanics such as aerodynamics, 
f-^ ^ combustion, nuclear fusion reactors or bioreactors. The presence of this phenomenon in these devices may be beneficial 
I— I , or detrimental, depending on each particular application [ll|, [l^ ■ 

Vortex breakdown (VB) has been firstly observed in open flows [l|, but subsequently it has been obtained also in 
' I experiments performed in confined fiows, for example, in closed cylinders 0]. It is worth noting that the characteristics 
• of VB in both cases are strongly similar, suggesting the possibility that the basic mechanism of VB is the same in 
I ■ both situations [To| . Experimental measurements show that flows in open channels can be accurately described with 
fS| • the relatively simple q- vortex model [5]. On the other hand, an analytical description for closed flows which are 
' considerable more complex than closed flows is not known [1, @, Q- Thus, experiments performed with flows confined 
, in cylinders are very attractive in order to understand the basic mechanism of VB, as well as investigate possible 
^-H ' control mechanisms. 

00 , Due to the great number of practical implications of VB, the development of mechanisms for controlling its emer- 
' gence is of considerable interest. Recently, different methods for controlling VB in closed fiows have been proposed 
^ ' using different techniques. The effect of adding a small, rotating or stationary, cylinder was first investigated using 
numerical and experimental techniques, by MuUin et al 1^. Husain et al [ll] considered the addition of co- or counter 



rotation near the axis using a small central rod rotating independently of the bottom or rotating end wall. They 
^ : concluded that such method may be an effective way to either enhance or suppress VB. Other app roaches are based 
■ - - ■ on the introduction of a small_ rotating disk in the end wall opposite to the rotating wall [l^ Il9l. [20j. or the use of 
axial temperature gradients [21|. 

In the Ref. ,12], the authors found that a straight small cylinder does not change qualitatively the characteristics 
of the fiow. They also considered the case of a stationary cylinder for fixed aspect ratio and Reynolds number they 
did not find important qualitative changes with respect to the case without inner cylinder. They concluded that it 
is necessary to slope the inner boundary in order to significantly alter the VB bubble. However, as these authors 
recognize it is difficult to find a reliable measure for the onset of VB using dye visualization. It is clear that a more 
accurate technique is necessary in order to get quantitative experimental results for the VB onset. 

The aim of this work is to analyze accurately the effect of the presence of a rod at the cylinder axis on the VB for an 
important range of parameters. Since the onset of VB is associated with the appearance of a stagnation point near the 
axis an accurate determination of the velocity field is necessary. Using the Digital Particle Image Velocimetry (DPIV) 
technique we obtain the critical Reynolds number for the onset of VB as a function of the aspect ratio and radius 
of the inner cylinder. In order to gain a deeper understanding of the internal structure of the flow and accuracy in 



2 



Rotating disc 







2d 


r— 




z 


<■ > 

R 



FIG. 1: The experimental setup consists of a closed cylinder, radius R and height H, with a rotating top wall. Along the axis 
of the cylinder a fixed rod of radius d is located. A external container was used in order to minimze the optical distortion of 
the images due to the curvature of th e cylindrical wall. 

the critical Reynolds numbers, we considered necessary to complement our experimental observations with numerical 
simulations. It is found that the existence of rods can affect significantly the critical value of the Reynolds number, 
Rc for the onset of the VB . 

This paper is organized as follows. In Sec.[TT]we present the experimental setup and general observations. In section 
III, we describe the numerical method used in the simulations. The comparison between the experiments and the 
numerical calculations is given in Sec. IIVI Finally, in Sec.|Vl we present a summary and the conclusions. 

II. EXPERIMENTAL SETUP 

The experimental setup consists of an plexiglas cylindrical container of inner radius R = 40.0 mm, a rotating disk 
at the top and a fixed disk at the bottom of the cylinder (see figure [J). The rotating disk is driven by an electronically 
controlled DC motor, operating between 17 = 1 rad/s and = 21 rad/s, with an error of less than 0.5% in the 
measurment of fl. The height of t he flow domain, and therefore the aspect ratio H / R, can be varied by moving 
the rotating disk up and down to a predetermined location. In this work, four different aspect ratios H / R: 1, 1.5, 2 
and 2.5 were used. 

In a first series of experiments, the emergence of vortex breakdown in absence of rod for comparison with results of 
previous works is studied. In the second part of the experimental work, the effect of axial rods on VB development is 
investigated. We use the cylinder described above, with an axial rod fixed on the bottom disk. The rod was perfectly 
aligned along the axis of the cylinder and was kept at rest during the experience. Three different rod sizes, with 
corresponding radius d = 1.0 mm, d = 2.5 mm and d = 5.0 mm were used. 

The fluid used was water dissolutions of glycerin at 60% in mass, with v = 1.0 x 10^^ m^/s. Neutrally buoyant 
polyamide particles of 50 x 10~^ m diameter was seeding inside fluid in o rder to obtain the velocity profile. A vertical 
green laser sheet of 100 mW, with 2 mm of width was used as an illumination system. The sheet was focalized at 
the meridional plane of the cylinder. The images of the flow were record using a CCD camera at 90 fps. In order to 
minimize optical distortion of images due to the curvature of the cylindriln the onset of VB a stagnation point in the 
axis appears in the flow . In this case, the vertical velocity at the axi s is null and we used this fact to 

Figure [2] shows the typical DPIV images for the onset of VB corresponding to different radius of the rod.cal wall, 
the container was immersed in a rectangular box filled with the same working fluid, since both the solution and the 
plexiglas have similar refractive indices. The experiments were conducted in an air-conditioned room at 20°C and 
the fluid temperature was monitored regularly using a thermocouple located at the bottom of the external container. 
The Reynolds number was deflned as in previous studies as Re = VlR^ /v, which varies in the experiments between 
600 and 2900, with a 2% of error. 

The conventional Digital Particle Velocimetry Velocity (DPIV) technique was implemented in order to obtain the 
velocity fleld inside the flow. This technique allow us to get accurate quantitative values of the velocity proflle based 
on the cross-correlation of two consecutive images. I n fact, this technique allow us deteminate accurately the critical 
Reynolds number (Rec) at the onset of VB through the measurment of the vertical velocity at the axis of the cylinder. 
The onset of VB is typically associated with localized regions of reverse flows, i.e. an internal stagnation point on the 
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FIG. 2: Magnifications of tlie zone corresponding to tfie recircuiation bubbie at tlie onset of VB for different radius of tfie rod. 
Tlie aspect ratio {H/R = 2) is tlie same in all the images while the Rcc varies in each image: (a) Re — 1474, without rod. (b) 
Re = 1313, d = 1.0mm, (c) Re = 1431, d = 2.5mm. (d)Re = 1671, d = 5.0mm. 
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FIG. 3: Axial velocity profiles at the axis (r = 0) for H/R — 2 obtained using the DPIV technique and corresponding to the 
cylinder without rod . (a) situation before the onset of VB, Re = 1391, the axial velocity at the axis is allways positive, (b) 
first VB, Re = 1474, the axial velocity at the axis present a null value at the stagnation point located at Zc = 3.1 cm. (c) onset 
of the second VB, Re — 1836, in this case the axial velocity at the axis presents a range of negative values in the place of the 
first reciculation bubble and a second stagnation point at Zc — 4.5 cm where the second bubble is going to appear. 



vortex axis appears, located at a certain height, (zc). The appearance of VB bubbles is a manifestation of internal 
separation of the streamlines from the axis, causing the reversal of the axial velocity. We used the criteria that VB 
arises when the axial velocity changes the sign at the axis. In the situation when a fixed rod is localizate at the axis 
of the cylinder, we analyze the vertical velocity at the wall of the rod. 

Figure [2] shows the typical DPIV images for the onset of VB corresponding to different rod radius. 

In figure [3] three typical axial velocity profiles at the axis (r = 0) for H/R = 2 correspondig to the cylinder whitout 
rod are show. In Fig. [3] (a) previous to the onset of VB the axial velocity is always positive and in this case the 
Reynold number is Re = 1391. At the onset of VB, Fig.[3](b) corresponding to Rcc = 1474, the axial velocity present 
a null value at the stagnation point located at Zc = 3.1 cm above the bottom of the cylinder. Finally, in Fig.[3](c) we 
can observe that the second VB take place at Rcc — 1836. In this case the axial velocity at the axis presents a range 
of negative values in the place of the first reciculation bubble and a second stagnation point at Zc = 4.5 cm above the 
bottom of the cyhnder. 
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III. NUMERICAL METHOD 



The numerical simulations considered here were obtained with the in-housc flow solver caffaSd.MB developed jointly 
by Universitat Rovira i Virgili (Tarragona, Spain) and Universidad de la Repiiblica (Montevideo, Uruguay). It is an 
original Fortran95 implementation of a fully implicit finite volume method for solving the 3D incompressible Navier- 
Stokes equations in complex geometry, using block structured body fitted grids. This three-dimensional solver is based 
on a two-dimensional solver described in |23| • Further description of this model can be found in Ref. [1^, [2^ . 




The mathematical model comprises the mass ([T|) and momentum ([2]) balance equations for an incompressible 
Newtonian fluid: 



These equations hold in any portion of the domain, being S the boundary of f2 and fis the outward normal vector 
at the boundary S. The momentum balance equation ^ has been expressed for the flrst component u of the velocity 
vector V = {u,v,w), with similar expressions holding for the other components. The viscosity /i of the fluid and the 
symmetric deformation tensor D were used for the viscous term. 

Since the Reynolds number was relatively low for all cases, since Re — 800 to Re = 2700, no turbulence model was 
required so transient solutions were computed directly. The time step was set to 10^^ s for all cases. Simulations 
were run starting from null velocity fields through 10^ time steps, or about 100 s of flow time. 

Two different grids were used, one to study the fluid without rods, and other to study the system including the 
rods. To study the system with rod C-grid block, with a no-shp boundary condition at walls. Fig. [D^b), is used. The 
different choice of grid topology for the two cases reflects that a boundary condition near the axis exists for the case 
that includes rods, while no boundary condition exists at the axis for the case without rods. Thus, the two-block 
grid strategy adopted for the second case ensures that the critical region near the axis is treated entirely as an inner 
region in the absence of a rod, not requiring the specification of delicate boundary conditions there. On the other 
hand, for the case that includes rods, the single cylindrical block is the natural choice allowing a direct specification 
of the non- slip boundary condition at the rod. On the other hand, the grid used for the cylinder without rods is a 
block structured grid, made of two blocks. The inner one is a prism of a square base of 5 mm (0.2592 x 10^ cells) 
located in the axis of the cylinder, while the outer one is a C-grid block that wraps around the first one, see Fig. [3] 
(a) (2.304 X lO^cells). 

The numerical model allows us to get the three components of the velocity field (w^, wg, Vz) and its derivatives. 
As mentioned before, the reversals of the axial component of the velocity Vz in the region near the axis gives us 
information on the emergence of areas of recirculation (bubbles), which are a way of determining VB. In Fig. [5] 
numerical results of the axial component of the velocity are showed for the cylinder without rod. The two showed 
panels correspond to one (top) and two (bottom) VB bubbles. This velocity profiles show an excellent agreement 
with those obtained experimentally (Fig. [3]). 

The simulations are performed for different cases (three values of the aspect ratio H/R, cylinder with and without 
rods) and Reynolds numbers. The numerical results are in excellent agreement with the experimental velocity profiles. 
The position of the stagnation point is obtained with a difference less than 5% in the system with and without rods. 
The numerical results were analyzed trough streamlines tp and azimuthal vorticity rye = rjy cos O — rj^ sin 6. In addition, 
the energy density, defined as e = ^ -I- ^(w^ + Vg + v^), the velocity module and the pressure field were obtained (not 
shown in this work). 

Figure [6] shows contours for the azimuthal vorticity rye in the meridional plane for H/ R — 2 a,t critical Reynolds 
number Rcc as indicated. The contour levels are non-uniformly spaced, with 20 levels in the positive range and 20 in 
the negative one determined by levelpos{i) — Max{variable) x {^)^ and levelnegi}) = Miniyariahle) x (^)^ @- All 
are plotted &i t — 100 s, when steady-state fiow conditions have been reached. 
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FIG. 4: (a) Cross section of the block-structured grid, with a central square block, used in the case of the cylinder without 
rod. This grid topology prevents the need for internal boundary conditions at the axis of the domain for the case without rod. 
(b) Cross section of the structured grid for the case of the cylinder with rod. 
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FIG. 5: Numerical results of the axial velocity profiles (left) and streamlines (right) corresponding to the same situation of 
figure [3] (a) and (b) one bubble is observed for H/R = 2 and Re = 1474. (c) and (d) two bubbles are observed for H/R — 2 
and Re = 1836. 



IV. RESULTS AND DISCUSSION 

For convenient comparison with previous works 0] , we first studied the emergence of vortex breakdown in absence 
of rods. In table llVl the experimental and numerical results for the onsent of VB for different aspect ratios are shown. 
In this table also Escudier results are shown. First, if we compare our experimental and numerical results, we can 
observe that its present an excellent agreement with a difference lower than 2%. In second place, our results are 
in very good agreement with the results of Escudier Q. Our results present a difference lower than 3% respect to 
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FIG. 6: Numerical results for the azimuthal component of the vorticity at the onset of VB and aspect ratio H/R = 2. The 
lined regions represent the fixed rod. (a) without rod, Re = 1474, (b) fixed rod of radius d = 1.0 mm and Re = 1340. (c) fixed 
rod of radius d = 2.5 mm and Re — 1435. (d) fixed rod of radius d = 5.0 mm and Re = 1700. The scale, indicated in the 
colorbar at the right, is the same for all the panels. 



Escudier's results. We can conclude that the method experimental and mumerical that we have used are very reliable. 
Thus, the DPIV technique provide an accuracy method to determinate the onset of the VB through the analize of 
vertical velocity profile. The result presents in this work shown that the DPIV technique are more realible that one 
based in visualization using dye technique. 

Now, we consider the situation in which an axial fixed rod of radius d is present. In table IIVI for d/R = 0.025 we 
can obaserve that the Reynolds number necessary to achive the onset of VB is less than the case without fixed rod. 
Meanwhile for d/R = 0.0625 the onset is achived for Rec slightly lower than d/R — and for d/R = 0.1250 the Rec 
present values 15% above than the case without rod. 

In figure Fig [7] is depicted a comparison between the experimental and numerical results concerning the effect of the 
cylinders on the first bubble formation, showing that there is a very good agreement. We show the Rcc as function 
d/R for the three aspect ratio that we have studied. In order to analyze the response of the system when the radius 
of rod is increased, we calculated numerically the Rcc for d/R — 0.25. The difference in the Reynolds number at the 
onset is approximately 65% between cylinder without and with rod, for all aspect ratio analyzed here. In summary, 
we can conclude that if the aspect ratio is kept constant, increased the radius of the rod represent a significative 
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TABLE L Experimental and numerical critical Reynolds numbers corresponding to the appearance of the first VB for different 
aspect ratios and radius of the rod. The experimental values (1) were taken from |3| and the (2) were obtained in this work 
using the DPIV technique. 
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FIG. 7: Experimental (open symbols) and numerical (filled symbols) values of the critical Reynolds numbers corresponding to 
the appearance of the first bubble as a function of the rod radius, for different values of the aspect ratio H/R = 1.5 (diamonds), 
H/R = 2.0 (circles) and H/R = 2.5 (squares) 



change in the Rcc at the onset and for small radius of the rod, the Rcc is lower than the case without rod. However, 
there is a region of radius of the rod, that depend of aspect ratio, which slightly variation is observed for the Rcc 
between cylinder with and without rod. 

A clear modification is noted also in the vorticity structure, as can be seen from Fig.[6l For d — 2.5 mm and d = 5.0 
mm, we observe that the critical Reynolds number to begin vortex breakdown is increased. The effect in the structure 
of the vorticity in these cases are similar to the observed for the smallest rod. Experimental and numerical for the 
critical value Rbc are summarized in Table [TVl 

It is interesting to compare our results with those obtained in references p^ . pdj . In these two works, it was 
concluded that a very slender rod at rest does not introduce significant changes in VB phenomena. The radius 
ratio used by these authors was d/R = 0.04, which is similar to the small radius ratio value we considered, that is 
d/R — 0.025. In reference [l^ . the authors used an inner cylinder with ratio d/R — 0.1. In our experiments we 
observed appreciable changes in the critical value of Re, for all values of the ratios considered. Looking carefully at 
Figures 2 and 3 of [l^, we observe slight differences which, according to the authors, are within the experimental 
errors. For this reason, we took specially careful measurements in these cases. 
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At this point it is interesting to consider the origin of the effect of the rod on VB, and in particular if it is of viscous 
nature or not. Inviscid theory has largely used to explain the origin of VB. The results has been sufficiently succesfuU 
to support the belief that VB may occur with no presence of viscosity. In ref . 7] , Brown and Lopez argumented that 
VB in closed flows takes place in conditions for which the flow may be considered as inviscid. In this case, a change 
of sign in the azimuthal component of vorticity rjg is expected to occurs near the stagnation point |7|] • This condition 
is verified in the VB that takes place without rod, as can be seen in figure [6^ • In this figure, the stagnation point 
occur at z = 0.032, where z measures the height from below. However, this condition is not satisfied in the presence 
of the rods. In these cases, the stagnation point occur near z = 0.03 on rod surface. This point is surrounded with 
a region of constant rjg sign, as can be seen in figures [6}D,c,d. Thus the relation of Brown and Lopez does not holds, 
showing that we cannot expect to describe the VB appropriately with only inviscid theory in the presence of rods. 

The decrease of Rcc due to small rods appears to be of viscous origin. The presence of the rods tends to reduce 
the velocity of the fluid due to stress at the rod surface, favoring the appearance of the stagnation point and thus 
reducing the critical value Rcc needed for the appearance of VB. The reducing effect of the viscosity is also present 
in the absence of rods and may begin the process that leads to the formation of the stagnation point (as occur in 
pipe flows f?\), but the presence of the rigid limit imposed by the rod enhances this effect very much. Thus, a very 
thin rod produces an significative effect, because the boundary conditions at the axis changes drastically from free to 
no-slip. This produces an appreciable change in Rcc. A similar result was obtained in previous works where it has 
been shown that even though a thin wire may promote the formation VB in delta wings [22]. On the other hand, we 
note that when the diameter d is sufficiently large, the critical value of Re increases with increasing d (see Fig. [7]). 
This effect is probably of inviscid origin. We support this conjecture with the results of a simple model for flows in a 
pipe, described in the appendix. With this model we obtain that the suppressing effect of increasing d also appears in 
unconfined pipe flows. Then, when d increases, this effect of inviscid nature compensates the viscous one and finally 
the VB is displaced to higher values of Re. 

V. SUMMARY AND CONCLUSIONS 

In this work, we studied the changes in the VB phenomena that are produced when a rod is located at the axis of 
the cylinder container i.e. the vortex centcrline. The experiments shown that the rod may increase or decrease the 
critical Reynolds number for the emergence of VB. The effect of small rods is to decrease Rcc, while Rc is increased 
when d/R is beyond a certain critical value [d/ R)c, which is of order 0.06. This critical value depends weakly on 
the aspect ratio, decreasing smoothly with increasing H/R. The effects of the rods of practical interest since may be 
used to perform a control of the VB in a way that is simpler than other previously proposed in the literature, since 
it does not require additional auxiliary devices. It is worth noting that this method, unlike the approach proposed in 
previous works, does not imply the addition of swirl near the axis. This simplicity is an interesting feature, making it 
more feasible to be used in engineering devices. The modification of the critical Re are very noticeable. The volume 
ratio (cylinder to row) is V1/V2 ~ 4 x 10"'^, while the decrease of the critical Reynolds number is about 10% . So 
that the shift of the critical Reynolds number is 20 times larger that the percent modification of the volume of the 
cylinder, showing the effectiveness of the method. The dependence of the critical Re with the rod diameter obtained 
experimentally is in very good agreement with the found with numerical calculations. Using a simple theoretical 
model for flow in pipes, we obtained that the influence of the inner cylinder on the VB works in similar fashion in 
open flows, mainly for not too small values of {d/R)c- We gave arguments to support that the enhancing effect of 
very small rods on VB is of viscous origin. 

We acknowledge financial support from the Programa de DesarroUo de Ciencias Basicas (PEDECIBA, Uruguay) 
and Grants FCE 9028 and PDT54/037 (Conicyt, Uruguay). 

APPENDIX 

In this appendix, we consider a model for VB in open swirling pipe flows. Our interest is to study how the VB 
is affected when an inner cylinder is located at the axis of the duct. With this purpose, we use a simple model 
based on the failure of the cylindrical solution. In the past, it has been shown that simple theoretical vortex models 
of open flows may predict with some accuracy the e merg ence of vortex breakdown. A classical approach based on 
axisymmetric inviscid analysis was given by Batchelor [23| . In this model, it is assumed that initially the fluid moves 
inside a cylindrical duct of radius ai , with a Rankine vortex velocity field, and goes to another duct limited by 
two cylindrical boundaries of different radius. When the flow develops the stagnation point, the radial component 
of the velocity becomes important and the cylindrical approximation is no longer valid [i*]. This is reflected in the 
disappearance of cylindrical solutions. Thus in this model, the disappearance of cylindrical solutions is identified with 
the emergence of VB [2^ . 
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In a cylindrical system of coordinates (r, 0, z) the velocity field of the assumed upstream flow is written as: 

Vr = 0, (3) 

\(jr < r < a2, , . 

= <„, ^ ^ (4) 
yV /r a2 < r < ai, 

= Uo. (5) 

where a,T and [/q are constants, with T — <ja\, in order to impose the continuity of the velocity. We introduce the 
streamfunction ^ defined as 

Vz = -^: Vr = (6) 

r or r oz 

Employing the momentum conservation we obtain the following equation for the streamfunction: 

9,1 a*, ^de ^dK 

r — f )H = r K (7) 

dr^r dr' dz^ rf* d* ^ ' 

where e = + ^ and K — rvg. For the flow given by Eqs. ([3]), we have that 

20-2 1 2ct 
e = — * + -Ul K = — *, /or < r < ai (8) 

The same relations hold downstream for the steady flow, and so that Eq. ([7]) takes the form for the rotational part, 
if we restrict ourselves to consider cylindrical solutions, 



*(r) = -t/or2 + AFi(7r) + BY^{-ir) (9) 

where Fi and Yi are the Bessel functions of the first and second kind respectively and ^ — 2(t /Uq [1^. The irrotational 
part of the flow is given hy vq — T /r, Vz — U, Vr — 0. 

We particularize now to consider the situation in which the upstream flow goes two a duct limited by to cylindrical 
surfaces of inner radius d and outer radius R. The unknown constants which appear in Eq. ^ can be determined 
imposing the mass conservation and the continuity of the pressure. From these conditions, we obtain the following 
equation that gives the value of the radius of the rotational core in the downstream duct c 



with 



^"'""^^ ^''lMlb)-^jYoi,c)-l = 0, (10) 



^ Up di4^c^)Y,i^d) + cd^Y,ijb) 
" 2cd^Ji(7c)Fi(7rf)- Ji(7rf)yi(7c)^' ^ ^ 



_ Up dial-c^)Mjd) + cd^Mjb) 
" 2cd^YiMMjd)-Y,ijd)M^cy- ^ ^ 

In figures [5] and [5] it is shown the value of c as a function of the nondimensional swirl parameter S, defined as 
S — "fR. From these figures we can see that at certain critical values Sc of S, two branches of solutions collide and 
disappear. According to the failure of the cylindrical approximation criteria, this is the value of S for which VB takes 
place. From these figures, we can see that the effect of a very slender inner cylinder may be to increase or diminish 
slightly the critical value Sc (the first case is observed in figure [5] and the second in figure [5]) . However, for all the 
situations considered, Sc increases with d, as long d is above a threshold. In this case, the emergence of VB was 
transferred to higher values of the nondimensional swirl parameter S and the VB was suppressed in the range of S 
values contained between the old and new critical values of Sc- These results shows that the present inviscid model 
predicts the suppressing effect of the slender cylinders, in analogy with the experimental observations of the closed 
flow, showing that viscosity is not necessary to explain this behavior. However, the decreasing effect on Sc that is 
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FIG. 8: Dimensionless core size, c/ R, as a function of the swirl, S, for different radius of the fixed rod, (a) d = 0, (b) d = 0.2, 
(c) d = 0.4, (d) d — 0.6 and (e) , d = 0.8. Other parameter values: Ro — 1, co — 0.4 and R — 1.2 



c/R 




10,0 



FIG. 9: As figure Ebut with Rq = 1, cq = 0.25 and R = 1.2, for (a) d = 0, (b) d = 0.2, (c) d = 0.4, (d) d = 0.6. 



observed for small values of d in figure [9l has nothing to do with the observed in the closed flow for the smallest of 
value of d, where the effect appears to be of viscous origin. 
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